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ABSTRACT 

We use the combination of the 2 Ms Chandra X-ray image, new J and H band images, and the Spitzer 
IRAC and MIPS images of the Chandra Deep Field-North to obtain high spectroscopic and photometric 
redshift completeness of high and intermediate X-ray luminosity sources in the redshift interval z = 2 — 3. 
We measure the number densities of z = 2 — 3 active galactic nuclei (AGNs) and broad-line AGNs in 
the rest- frame 2 — 8 keV luminosity intervals 10^"' — 10'*^ and 10'*'^ — 10'*'' ergs s^^ and compare with 
previous lower redshift results. We confirm a decline in the number densities of intermediate- luminosity 
sources at z > 1. We also measure the number density of z = 2 — 3 AGNs in the luminosity interval 
lO''"^ — lO'*'''^ and compare with previous low and high-redshift results. Again, we find a decline in the 
number densities at z > 1. In both cases, we can rule out the hypothesis that the number densities 
remain flat to z = 2 — 3 at above the 5cr level. 

Subject headings: cosmology: observations — galaxies: active — galaxies: evolution — galaxies: 
formation — galaxies: distances and redshifts 
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1. INTRODUCTION 

Low-redshift hard X-ray luminosity functions have been 
well determined from the combination of highly spectro- 
scopically complete deep and wide-area Chandra X-ray 
surveys. At z < 1.2, the hard X-ray luminosity functions 
for active galactic nuclei (AGNs) of all spectral types and 
for broad-line AGNs alone are both well described by pure 
luminosity evolution, with evolving as (1 -I- z)"^'^ and 
(1 + z)^ °, respectively (Barger et al. 2005). AGNs decline 
in luminosity by almost an order of magnitude over this 
redshift range. 

Barger et al. (2005) compared directly their broad-line 
AGN hard X-ray luminosity functions with the optical 
QSO luminosity functions from Groom et al. (2004) and 
found that the bright end luminosity functions agree ex- 
tremely well at all redshifts. However, the optical QSO 
luminosity functions do not probe faint enough to see 
the downturn in the broad-line AGN hard X-ray lumi- 
nosity functions at low luminosities and may be missing 
some sources at the very lowest luminosities to which they 
probe. 

The Groom et al. (2004) pure luminosity evolution is 
slightly steeper than that of Barger et al. (2005), but 
within the uncertainties, the two determinations are con- 
sistent over the z — — 1.2 redshift interval. To investi- 
gate whether pure luminosity evolution continues to hold 
at higher redshifts, Barger et al. (2005) used the Groom et 
al. (2004) evolution law (which was fitted over the wider 
redshift range z = 0.3 — 2.1) to correct all of their broad- 
line AGN hard X-ray luminosities to the values they would 
have at z = 1. They then computed the broad- line AGN 
hard X-ray luminosity functions over the wide redshift in- 
tervals z = 0.2 - 0.7, 0.7 - 1.5, and 1.5 - 2.5. Barger et 
al. (2005) found that the lower redshift luminosity func- 
tions matched each other throughout the luminosity range. 



while the highest redshift luminosity function matched the 
lower redshift functions only at the bright end, where the 
optical QSO determinations were made. They therefore 
concluded that there are fewer intermediate X-ray lumi- 
nosity broad-line AGNs in the z = 1.5 — 2.5 redshift inter- 
val, and hence that the pure luminosity evolution model 
cannot be carried reliably to the higher redshifts. This 
was consistent with other analyses made of less complete 
samples of all spectral types together, which had found 
evidence for peaks and subsequent declines in the num- 
ber densities of both intermediate-luminosity (e.g., Cowie 
et al. 2003; Barger et al. 2003a; Hasinger 2003; Fiore et 
al. 2003; Ueda et al. 2003) and high-luminosity (Silver- 
man et al. 2005) sources, with the intermediate-luminosity 
sources peaking at lower redshifts than the high-luminosity 
sources. 

Recently, Nandraet al. (2005; hereafter, N05) have ques- 
tioned the evidence for a decline in the number densities 
of intermediate-luminosity AGNs at z > 1. By combin- 
ing deep X-ray data with the Lyman break galaxy (LBG) 
surveys of Steidel et al. (2003), they argue that the num- 
ber densities of intermediate-luminosity AGNs are roughly 
constant with redshift above z = 1. 

Only with extremely deep X-ray data and highly com- 
plete redshift identifications can we address the above con- 
troversy about a decline. In this paper, we measure the 
z = 2 — 3 number densities of both high and intermediate 
X-ray luminosity sources in the 2 Ms Chandra Deep Field- 
North (CDF-N). In addition to our highly complete spec- 
troscopic redshift identifications, we also use the recently 
released Spitzer Great Observatories Origins Deep Survey- 
North IRAC and MIPS data, in combination with deep J 
and H band data obtained from the ULBCAM instru- 
ment on the University of Hawaii's 2.2 m telescope, to es- 
timate accurate infrared (IR) photometric redshifts for the 



^Department of Astronomy, University of Wisconsin-Madison, 475 North Charter Street, Madison, WI 53706 
^Department of Physics and Astronomy, University of Hawaii, 2505 Correa Road, Honolulu, HI 96822 
^Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

1 



2 



X-ray sources that could not be spectroscopically or op- 
tically photometrically identified. We assume = 0.3, 
flA = 0.7, and Hq = 70 km s~^ Mpc"^ throughout. 



2. DATA 

The deepest X-ray image is the w 2 Ms CDF-N ex- 
posure centered on the Hubble Deep Field-North taken 
with the ACIS-I camera on Chandra. Alexander et al. 
(2003) merged samples detected in seven X-ray bands 
to form a catalog of 503 point sources over an area of 
460 arcmin^. Near the aim point, the limiting fluxes 
are « 1.5 x 10~^^ ergs cm~^ s~^ (0.5 — 2 keV) and 
w 1.4 X 10"^^ ergs cm'^ s"^ (2-8 keV). 

The ground-based, wide-field optical data of the CDF-N 
are summarized in Capak et al. (2004), and the HST ACS 
GOODS-N data are detailed in Giavahsco et al. (2004). 
The Spitzer IRAC and MIPS data are from the Spitzer 
Legacy data products release and are presented in Dick- 
inson et al. (2005). The new ULBCAM J and H band 
data, which cover the whole CDF-N area to la depths of 
25.2 in J and 24.2 in H (these are measured in 3" aper- 
tures and corrected to total magnitudes), are described in 
Trouille et al. (2005). All magnitudes in this paper are 
in the AB system. Most of the spectroscopic rcdshifts 
for the X-ray sources are from Barger et al. (2002, 2003b, 
2005), but we also obtained some additional rcdshifts with 
the DEIMOS spectrograph on the Keck II telescope dur- 
ing the Spring 2005 observing season. One spectroscopic 
redshift {z = 2.578) is taken from Chapman et al. (2005). 



2.1. X-ray Incompleteness 

We consider in our analysis two restricted, uniform, fiux- 
limited X-ray samples. The most important consideration 
in choosing our soft X-ray flux limits was that there be 
no signiflcant X-ray incompleteness. In Figure 1, we show 
0.5 — 2 keV flux versus off-axis radius for the full 2 Ms 
CDF-N X-ray sample of Alexander et al. (2003) {squares 
and diamonds). We use small squares to denote sources 
that are detected in the 2 Ms survey but not in the 1 Ms 
survey (Brandt et al. 2001). For the bright subsample, we 
consider sources in a 10' radius region with fluxes above 
1.7 X 10"^^ ergs cm'^ s"^ (0.5 - 2 keV), and for the deep 
subsample, we consider sources in an 8' radius region with 
fluxes above 1.7 x lO"^*^ ergs cm'^ s~^ (0.5 - 2 keV). We 
show these flux and radius limits with solid lines. The two 
subsamples comprise 160 sources, 8 of which are stars. 

We can see from Figure 1 that the flux limits of our sub- 
samples arc well above the flux limit of the 2 Ms survey 
{dotted curve) at all radii, suggesting that there should be 
very little X-ray incompleteness in our sample. Indeed, 
when we consider the 160 sources in the 2 Ms exposure 
that lie within our two subsamples, we flnd that only five 
sources were not already detected in the 1 Ms exposure. 
If we were using the 1 Ms catalog in our analysis rather 
than the 2 Ms catalog, then this would correspond to a 
3% incompleteness. With the 2 Ms exposure, we expect 
X-ray incompleteness to be negligible. 
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Fig. 1. — 0.5 — 2 kcV flux vs. off-axis radius for the fuU 2 Ms 
CDF-N sample of Alexander et al. (2003) {squares and diamonds). 
The flux and radius limits of our bright and faint subsamples are 
denoted by the solid lines, and the 2 Ms flux limit is denoted by 
the dotted curve. Sources detected in the 2 Ms survey but not in 
the 1 Ms survey (Brandt et al. 2001) are denoted by small squares. 
Only five such sources lie within our subsamples. Sources wfith spec- 
troscopic redshifts in the interval 2 = 2 — 3 are denoted by large 
squares. Our flux and radius limits were chosen without reference 
to the distribution of the z = 2 — 3 sources. 

2.2. Redshift Interval 

Having selected a highly complete X-ray sample, our 
next concern was that we sample the full luminosity ranges 
of interest without clipping out any sources. For the tra- 
ditional intermediate and high X-ray luminosity cut-offs 
of L2-8 kcV = lO^'' and 10*^ ergs s^^, this very naturally 
sets the redshift interval to 2 = 2 — 3. We can see this 
from Figure 2, where we show rest-frame 2 — 8 keV lumi- 
nosity versus redshift for the CDF-N X-ray sources in our 
(a) bright and (b) faint subsamples. We determined the 
rest-frame 2 — 8 keV luminosities from the observed- frame 
0.5 — 2 keV fluxes, assuming an intrinsic F = 1.8 spec- 
trum. At 2; = 3, observed-frame 0.5 — 2 keV corresponds 
to rest-frame 2 — 8 kcV, providing the best possible match 
to lower redshift data. Moreover, the 0.5 — 2 keV Chandra 
images are deeper than the 2 — 8 keV images, so using the 
observed-frame soft X-ray fluxes at high redshifts results 
in increased sensitivity. We only show the luminosities at 
low redshifts for illustrative purposes, since normally one 
would use the observed-frame 2 — 8 keV fluxes to calcu- 
late the 2 — 8 keV luminosities at these redshifts. The soft 
X-ray flux limits of our faint and bright subsamples (solid 
curves ) correspond to rest-frame 2 — 8 keV luminosities of 
10*'' orgs s^^ and 10*'* ergs s^"'^, respectively, at z = 3. 
Thus, we will not be clipping out any sources in either the 
-t'2-8 keV = 10*^ to 10** ergs s~^ interval in the faint sub- 
sample, or the L2-8 kcV = 10** to 10*^ ergs s"""^ interval 
in the bright subsample, if we choose 2 = 2-3 for our 
analysis. 

Once we have chosen our soft X-ray flux and radius lim- 
its and our redshift interval, we can plot the spectroscop- 
ically identifled z = 2 — 3 sources in the flux-radius plane 
and see how many lie in our restricted subsamples. This is 
shown in Figure 1, where we denote with large squares the 
sources that have been spectroscopically identifled to lie 
in the redshift interval 2 = 2-3. The soft X-ray flux and 
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Fig. 2. — Rcst-framc 2 — 8 kcV luminosity versus rcdshift for the CDF-N X-ray sources in our (a) bright and (b) faint restricted subsarnplcs. 
The luminosities were determined using the observed-frame 0.5 — 2 keV X-ray fluxes and if-corrections calculated using a F = 1.8 spectrum. 
Sources with spectroscopic redshifts are denoted by small solid diamonds. Broad-line AGNs are denoted by large solid diamonds. Sources 
with IR photometric redshifts are denoted by open squares. The eight remaining sources without spectroscopic or photometric redshifts (one 
in the bright subsample and seven in the faint subsample) arc denoted by open circles and are shown at z = 2.5. The solid curve in (a) 
shows the bright subsample 0.5 — 2 keV flux limit of 1.7 X 10^^''' ergs cm^'^ s^-*^ in a 10' radius region, and the solid curve in (b) shows 
the faint subsample 0.5 — 2 keV flux limit of 1.7 X 10^"*^^ ergs cm^^ s^^ in an 8' radius region. The dotted horizontal line in (a) shows the 
high-luminosity interval lower limit of 10^^ ergs s~^, and the dotted horizontal lines in (b) show the intermediate-luminosity interval upper 
and lower limits of 10^"* and IC*^ ergs s~^, respectively. 



radius limits for our two subsamples were chosen without 
regard to this distribution, as can be seen from the figure. 

There are 15 sources in our two subsamples with spectro- 
scopic redshifts in the interval z — 2 — 'i. 

2.3. Probability of Misidentifications 

Another issue we need to address is how many of the 
X-ray sources with detected optical/near- infrared (NIR) 
counterparts we may be misidentifying. Of the 160 sources 
in the 2 Ms exposure that lie within our two subsamples, 
138 are detected above the 3(t limit of 23 in H. Thus, only 
22 of the sources in our subsamples are not detected at the 
3fT level in the NIR. By randomizing the source positions 
and remeasuring the H band magnitudes many times, we 
found an 11% probability of misidentification. Thus, of the 
138 sources in our subsamples with 3a H band detections, 
we would expect that about 15 might be contaminated by 
the projection of random, superposed objects. 

Of the 22 sources that are not detected at the 3(7 level in 
H , 11 are detected in the IRAC 3.6/xm band at < 23 mag- 
nitude, with a 4% probability of misidentification based 
on randomized measurements. This corresponds to about 
one-half of a source. Of the remaining 11 sources that are 
not detected at 3.6/im, 5 are detected at i? < 25.5, with a 
10% probability of misidentification based on randomized 
measurements. Again, this corresponds to about one-half 
of a source. 

Thus, in total, we have optical/NIR identifications for 
154 of the 160 sources, eight of which are stars, and we ex- 
pect about 16 (or just over 10%) of these to be false. Even 
if a large fraction of the 16 real sources suff'ering from ran- 
dom projections were in the z = 2 — 3 range, which is 
not very likely, the correction to the number densities we 
derive subsequently would be small. 

3. PHOTOMETRIC REDSHIFTS 



Barger et al. (2002, 2003b) computed photometric red- 
shifts for the CDF-N X-ray sources using broadband 

galaxy colors and the Bayesian code of Benitez (2000). 
They only used sources with probabilities for the photo- 
metric redshift of greater than 90%, resulting in about an 
85% success rate for photometric identifications. These 
redshifts are robust and surprisingly accurate (often to 
better than 8% when compared with the spectroscopic red- 
shifts) for non-broad-line AGNs. 

We now have the advantage of the Spitzer IRAC 3.6/xm, 
4.5/im, 5.8/im, and 8.0/im data, as well as the NIR J 
and H band data, for estimating photometric redshifts. 
Barger et al. (2005) classified the X-ray sources into four 
spectral classes: absorbers (i.e., no strong emission lines 
[EW([OII])< 3 A or EW(Ha-|-NII)< 10 A]), star formers 
(i.e., strong Balmer lines and no broad or high- ionization 
fines), high-excitation (HEX) sources (i.e., [NeV] or CIV 
lines or strong [OIII] [EW([OIII] 5007 A) > 3 EW(H/3)]), 
and broad-line ACNs (i.e., optical lines having FWHM line 
widths greater than 2000 km s^^). The measured spectral 
energy distributions (SEDs) of the X-ray sources in each 
of these spectral classes turn out to be remarkably tight. 
Thus, we can construct templates from the average SEDs 
for each class and obtain an IR photometric redshift by 
fitting the templates to the optical through mid-infrared 
(MIR) data. We note that the template fitting is insensi- 
tive to the spectral class for sources without strong AGN 
signatures. 

In Figure 3, we compare both (a) the optical photo- 
metric redshifts from Barger et al. (2003), and (b) our 
new IR photometric redshifts with the spectroscopic red- 
shifts. Both photometric redshift techniques tend to fail 
on some of the broad-line AGNs (large open squares), 
but, fortunately, broad- line AGNs are straightforward to 
identify spectroscopically, even in the so-called redshift 
"desert" at ^ ~ 1.5 - 2, and nearly all of the CDF-N X-ray 
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Fic;. 3. — (a) Comparison of the optical photometric redshift estimates with Bayesian probabilities > 0.90 from Baigci et al. (2003b) with 
the spectroscopic redshifts. (b) Comparison of the IR photometric redshift estimates from this paper with the spectroscopic redshifts. Large 
open squares denote broad-line AGNs. 



sources have now been spectroscopically observed. While 
the Bayesian optical photometric redshift technique gives 
slightly tighter values at low redshifts due to the presence 
of strong spectral features and breaks in the templates 
from Coleman, Wu, & Weedman (1988) and Kinney et 
al. (1996), it does not do so well on the z > 2 sources. 
In our subsequent analysis, we adopt the more robust IR 
photometric redshifts when a spectroscopic redshift is not 
available. In Figure 4, we show histograms of the spectro- 
scopic (light shading) and IR photometric (dark shading) 
redshift identifications versus R magnitude for the sources 
in our (a) bright and (b) faint subsamples. 

The photometric redshifts not only tell us which of the 
spectroscopically unidentified sources are likely to have 
redshifts within the z = 2 — 3 redshift interval, but they 
also enable us to remove contaminants that arc really at 
lower redshifts. In Figure 5, we show four example mea- 
sured X-ray source SEDs (solid squares) redshifted to the 
rest-frame using the printed spectroscopic or IR photo- 



metric redshift. We have superimposed on each SED the 
average spectral template chosen by the IR photometric 
redshift fitting for that source. We selected these four ex- 
amples to illustrate (a, b) two spectroscopically identified 
sources with different spectral classes (also printed), both 
of which are in the redshift interval z = 2 — 3, (c) one IR 
photometric redshift of a lower redshift source, and (d) one 
IR photometric redshift of a source in the redshift interval 
z = 2-3. 

Of the 44 sources with /0.5-2 keV > 1.7 x 
lO"""^^ ergs cm"^ s^^ within a 10' radius, 38 have spec- 
troscopic redshifts. Of the remaining 6 sources, 5 have 
a robust IR photometric redshift, none of which lie in 
the redshift interval z = 2 — 3. Of the 136 sources with 
/0.5-2 keV > 1.7x 10~^^ ergs cm~^ s~^ within an 8' radius, 
88 have spectroscopic redshifts. Of the remaining 48, 41 
have robust IR photometric redshifts, 12 of which lie in 
the redshift interval z = 2 — 3. That leaves 8 sources in 
total with neither a spectroscopic nor a photometric red- 




FiG. 4. — Fraction of the total (open histogram) (a) bright and (b) faint restricted subsamples that are spectroscopically (light shading) 
and photometrically (dark shading) identified vs. R magnitude. All R> 27 sources are placed in the final magnitude bin. 
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Fig. 5. — Bolometric flux (i^fi^) vs. rest frequency for four example X-ray sources in the CDF-N sample. Solid squares show the optical 
{HST ACS and ground-based), near-infrared (J and H), and mid-infrared (3.6/im, 4.5/im, 5.8/um, 8.0/im, and 24/im) measurements. Open 
squares denote 5a upper limits. Curves show the average spectral energy distributions of the CDF-N X-ray sources of a given spectral class 
selected as templates by the IR photometric redshift estimation code (see text for details). Spectroscopic redshifts and spectral classes are 
used in (a) and (b), and IR photometric redshifts are used in (c) and (d). 



shift; 6 were not clearly detected in the H, IRAC 3.6/zm, 
or R bands (see §2.3), and the photometric redshifts are 
not obvious from the SEDs for the remaining two. 

4. AGN NUMBER DENSITIES AT Z = 2 - 3 

Six of the eight sources without spectroscopic or photo- 
metric redshifts could, if they were assigned redshifts at 
the center of the redshift interval z = 2 — 3, have lumi- 
nosities that would place them into one of our two lumi- 
nosity intervals (one in the high-luminosity interval, and 
five in the intermediate- luminosity interval). In Figure 2, 
we denote the X-ray sources with spectroscopic redshifts 
by small, solid diamonds, and the broad-linc AGNs by 
large, solid diamonds. We use open squares to denote the 
sources with IR photometric redshifts, and we use open 
circles to denote the unidentified sources (these are plot- 
ted at 2: = 2.5). 

Wc computed the AGN number densities at z = 2 — 
3 for the spectroscopically or photometrically identified 
sources in the rest-frame 2 — 8 keV luminosity intervals 
1043 _ 1044 ergs s-i and 10^'' - 10"^^ ergs s"^ using the 
appropriate areas and volumes. These number densities 



are shown in Figure 6a as a solid diamond and a solid 
circle, respectively. The Poissonian la uncertainties are 
based on the number of sources in each redshift interval. 
Our high-luminosity sample consists of 8 sources in this 
redshift interval, all of which are broad-line AGNs. Our 
intermediate-luminosity sample consists of 13 sources in 
this redshift interval, none of which is a broad-line AGN. 

We also computed upper limits on the number densi- 
ties at z = 2 — 3 by assigning redshifts of 2.5 to the six 
unidentified sources described above. We denote these 
upper Hmits by solid (lO"*"^ — 10^^ ergs s~^) and dashed 
(1044 _ io45 gj-gg g-i) horizontal lines in Figure 6a. 

For comparison, we also plot in Figure 6a the AGN num- 
ber densities at z < 1.5 for our two luminosity intervals. 
These number densities were determined from the spec- 
troscopic sample of Barger et al. (2005), which includes 
the CDF-N (Barger et al. 2003), CDF-S (Szokoly et al. 
2004), and CLASXS (Stcffen et al. 2004) data. We do 
not expect cosmic variance to be an issue on these scale- 
lengths (Yang et al. 2005). We calculated the rest- frame 
2 — 8 keV luminosities for these sources using the observed- 
frame 2 — 8 keV band, assuming an intrinsic F = 1.8 spec- 
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trum. The Poissonian la uncertainties are again based 
on the number of sources in each redshift interval. The 
Barger et al. (2005) spectroscopic sample is substantially 
complete at these redshifts and luminosities, and we ex- 
pect that any incompleteness correction would be small. 

In Figure 6a, we see that the AGN number densities 
in both luminosity intervals show a steep rise at z < 1. 
The intermediate- luminosity number densities (diamonds) 
then show a marked decline to z = 2 — 3, while the 
high-luminosity number densities ( circles ) remain rela- 
tively constant to z = 2 — 3. If the number densities in 
the L2-8 kcV = 10^^ — 10^'' ergs s~^ interval were constant 
with redshift and equal to the peak value of ~ 10~^ Mpc~^ 
seen just below z = 1 in this luminosity interval, then we 
would expect 66 sources in the z = 2 — 3 interval at off- 
axis radii less than 8'. In fact, we only have 13 sources 
with spectroscopic and photometric redshifts in this red- 
shift interval in our faint subsample. Even if we were to 
add in all eight of the unidentified sources, this number 
would only rise to 21. Thus, we can reject the; hypothesis 
that the number densities do not decline with redshift at 
above the 5a level. 

In Figure 6b, we show the broad-line AGN number den- 
sities for the same two luminosity intervals. The high- 
luminosity number densities show a dramatic rise at z < 
1.5 (open triangles), and then a flattening to z = 2 — 3 
(solid triangle), while the intermediate-luminosity num- 
ber densities remain relatively constant at z < 1.5 (open 
squares), and then decline to z = 2 — 3, where we only have 
an upper limit (solid square and arrow). The z < 1.5 be- 
havior can be understood in the context of the broad-line 
AGN luminosity functions given in Figure 18 of Barger et 
al. (2005). Over the z = 0—1.5 redshift interval, the 
evolution of the broad-line AGN luminosity function is 



well defined by pure luminosity evolution with a rapid 
increase of luminosity with redshift. As a consequence, 
the 1/2-8 keV = 10**'* — 10**^ ergs s~* interval, which lies on 
the steeply declining high-luminosity end of the luminosity 
function throughout this redshift interval, rises rapidly. In 
contrast, the L2-8 kcV = 10^^ — 10'''* ergs s^^ interval lies 
just above the peak luminosity in the luminosity function 
at z = 0, and the number densities in this luminosity in- 
terval stay relatively constant over the z = — 1.5 redshift 
interval. 



5. COMPARISON WITH NANDRA ET AL. (2005) 

N05 combined deep X-ray data of the CDF-N (2 Ms) and 

of the Groth-Westphal Strip (GWS; 200 ks) with the LBG 
surveys of Steidel et al. (2003) to estimate the number 
density of intermediate-luminosity AGNs in the interval 
z = 2.5 — 3.5. Their method differs from the X-ray follow- 
up surveys in that it uses a rest-frame UV-selected sample 
over a narrow range of redshift with a cosmological volume 
determined from the optical selection function. The X-ray 
data are only used to determine whether an LBG hosts 
an AGN and to calculate an X-ray luminosity. The Groth 
Strip data are much shallower than the CDF-N data, and 
they apply an empirically determined 20% correction for 
X-ray incompleteness. 

To compare our results directly with those of N05, 
we used our CDF-N faint subsample to compute the 
AGN number densities at z = 2 — 3 in the same rest- 
frame 2 — 8 keV luminosity interval, L2-8 kcV = 10'*'' — 
-|^q44.5 gj.gg nsed by N05. We note that by choosing 
to use a broader luminosity interval that goes to higher 
luminosities, N05 is pushing into the QSO regime, where 
we know the number densities continue to rise to higher 
redshifts than the z = 1 redshift at which the intermediate- 
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Fig. 6. — (a) Number density of sources with rest-frame 2 — 8 keV luminosities between 10^'* and 10^^ ergs s^^ (diamonds) and 10^^ 
and IC*® ergs (circles). The solid symbols were determined using sources with both spectroscopic and IR photometric redshifts between 
z = 2 — 3, and the open diamonds and open circles at z < 1.5 were determined using the spectroscopic sample of Barger et al. (2005). 
Points below (above) z = 2 were determined from the observed-frame 2 — 8 keV (0.5 — 2 keV) sample. An intrinsic F = 1.8 was assumed, 
for which there is only a small differential Jf-correction to correct to rest-frame 2 — 8 keV. Poissonian la uncertainties are based on the 
number of sources in each redshift interval. The solid (dotted) horizontal bar shows the maximal number densities in the 10*^ — 10^* ergs 
{W^^ - 10*5 grgs s-i) interval found by assigning redshifts of 2.5 to the 5 (1) unidentified sources that could lie in these luminosity intervals 
at that redshift. (b) Number density of broad-line AGNs with rest-frame 2 — 8 keV luminosities between 10*^ and 10** ergs (squares) 
and 10** and 10* ergs (triangles). The symbols were all determined using sources with spectroscopic broad-line AGN classifications 
from Barger et al. (2005). The solid symbols were determined for the z = 2 — 3 interval. The solid square and horizontal bar show the la 
upper limit for the 10*^ — 10** ergs and z = 2 — 3 interval, since no broad-line AGNs were detected here. 
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luminosity number densities peak (see Fig. 6a). Our 
z — 2 — 3 data point (solid diamond in Fig. 7) contains 17 
sources, while the N05 data point contains 10. We com- 
puted an upper limit on our z = 2 — 3 number density by 
assigning rcdshifts of 2.5 to the five imidcntificd sources 
in our faint subsamplc that could lie in this luminosity in- 
terval if they were at that rcdshift. We denote this upper 
limit by a solid horizontal line in Figure 7. 

N05 claimed that their number density estimate in the 
z = 1.5 — 3 range was approximately an order of magni- 
tude higher than that of Cowie et al. (2003) and roughly 
equal to the Cowie et al. (2003) upper limit and a fac- 
tor of about 3 higher than that of Barger et al. (2005). 
N05's comparison with Cowie et al. (2003) was based on 
the L2-8 kcV > 10^^ ergs s^^ evolution, and hence was 
not shown on their Figure 1. They did show their esti- 
mates of the Barger et al. (2005) space densities, which 
they obtained by approximately fitting and then integrat- 
ing the Barger et al. (2005) luminosity functions over the 
kcV = 10^'' — 10^^"'' ergs range and then accoimt- 
ing for the differences in the adopted cosmology in the 
X-ray bandpass, but they did not show the upper limits 
from Barger et al. (2005) on their figure. 

For the purposes of a clear comparison between the 
different analyses, wc include on Figure 7 the Cowie et 
al. (2003) data points (solid squares) recalculated for the 
L2-8 keV = 10^^ — 10^^'^ ergs s~^ luminosity interval and 
the N05 cosmology (which is the same cosmology used 
throughout this paper). We also show the upper limit 
given by Cowie et al. (2003) for the z = 2 — 4 interval 
(dashed horizontal line). We denote the N05 result by 
a solid circle. We can see that there is no inconsistency 
between the Cowie et al. (2003) points and the present 
work, nor between the Cowie et al. (2003) upper limit and 
the N05 result. However, our present work has tightened 
up the measured number density and decreased the upper 
limit. 

We also computed the z < 1.5 AGN number densities in 
the 1/2-8 keV = 10^'^ — lO**^-^ ergs s~^ luminosity interval 
(open diamonds in Figure 7) using the spectroscopic sam- 
ple of Barger et al. (2005) . The decline in the AGN number 
densities at z > 1 is again highly significant. Since N05 
reached a diflferent conclusion, namely, that their "z = 3 
estimate is consistent with, and supportive of, the hypoth- 
esis that AGN activity remained roughly constant between 
z = 1 and z = 3" , we have also included on Figure 7 the 
Ueda et al. (2003) points (open triangles) and limits (dot- 
ted horizontal lines) that N05 compared with. The Ueda 
et al. (2003) values are generally similar to the Barger et 
al. (2005) values, except at z 1, where they have a high 
point. This is a consequence of the fact that the Ueda et 
al. (2003) data at this redshift come almost entirely from 
the CDF-N data, which are known to contain a redshift 
sheet at a median redshift of z = 0.94 (Barger et al. 2003). 
The Barger et al. (2005) number densities are smoothed 
out by the inclusion of the much wider CLASXS survey 
data (Steffen et al. 2004). 

We measure a value of 2.6to e >^ lO"'"^ Mpc'^ for the 
number density of sources in the L2-8 koV = 

10^3 - 

lO"''*-^ ergs s~^ range in the z = 2 — 3 interval, which may 
be compared with the N05 value of 4.2lJ;4 Mpc~^. In 
both cases, the errors only reflect Poissonian uncertainties 



from the small number of sources and are 68% confidence 
limits (Gehrels 1986). Although the N05 data point is 
statistically consistent with the present number density, 
the larger errors in N05 marginally permit a flat number 
density, while the present data deflnitively rule it out. 
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Fig. 7. — Number density of sources with rest- frame 2 — 8 keV lu- 
minosities between 10''^ and lO*'*'® ergs s~^. The sohd diamond at 
z = 2 — 3 was determined using sources with both spectroscopic and 
IR photometric redshifts, and the open diamonds at z < 1.5 were 
determined using the spectroscopic sample of Barger et al. (2005). 
The solid horizontal bar shows the maximal number density found 
by assigning redshifts of 2.5 to the 5 unidentified sources that could 
lie in this luminosity interval at that rcdshift. The solid circle at 
^ = 3 is from N05 from their LBG surveys. The solid squares at 
z = 2.5 and z = 3.5 arc from the spectroscopic sample of Cowie et 
al. (2003). The dashed horizontal line at 2: = 2 — 4 shows the max- 
imal number density for the Cowie ct al. (2003) sample, found by 
assigning rcdshifts to all of the unidentified sources in their sample 
that could lie in that rcdshift (and then luminosity) interval rcd- 
shifts at the center of the interval. The open triangles and dotted 
horizontal bars show the Ueda et al. (2003) points and upper limits, 
respectively. Poissonian la uncertainties are based on the number 
of sources in each redshift interval. 

While the present data unambiguously show that the 
number densities at the intermediate luminosities are de- 
clining to high redshifts, it should be noted that N05 were 
not correct in stating that a flat number density curve 
would argue against previous suggestions that the major- 
ity of black hole accretion occurs at low redshifts, around 
z = 1. The cosmological time available at high redshifts 
is much shorter than that at low redshifts, and the inte- 
gral of the energy density production rate over the cosmic 
time from z = to larger z would still be dominated by 
the sources around z = 1 for the case of a number density 
which rises to z = 1 and is constant at higher redshifts. 

It is worth briefly discussing the three major advantages 
that N05 argued their LBG method had over the X-ray 
follow-up method for determining AGN number densities. 
First, they claimed that the optical LBG selection func- 
tion is very well defined. However, as they noted in their 
paper, their volume element was calculated for galaxies, 
not AGNs, and the selection functions for both the broad- 
line AGNs and the narrow-line AGNs are likely to be quite 
different than that for the non-AGN LBGs. for which they 
referenced Steidel et al. (2002) and Hunt et al. (2004). 
They also noted that their method has the disadvantage 
of missing any AGNs that are too faint to be selected in 
their optical survey and/or have colors that fail the LBG 
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selection criterion. Hunt ct al. (2004) pointed out that 
three of the X-ray sources with spectroscopic redshifts in 
the range 2.5 < z < 3.5 reported by Barger et aL (2003b) 
were not picked up in the Steidel et al. (2003) LBG survey, 
and with the Barger et al. (2005) observations, that num- 
ber has now doubled to six. N05 argued that this incom- 
pleteness should be accounted for in their effective volume 
calculation (they only expect to pick up about 40% of the 
objects at 2 = 2.5 — 3.5 compared to a top hat function), 
but given that this is such a large correction, the uncer- 
tainties on the optical LBG selection function for AGNs 
are a concern. Direct selection of X-ray sources avoids this 
problem and should be much more robust. 

Second, N05 argued that they can apply a lower X-ray 
detection threshold for their subsample of LBGs than can 
be applied when constructing a purely X-ray based cata- 
log, thereby making their X-ray detections more complete. 
They declared this to be critical, since they are dealing 
with a population of sources which are close to the de- 
tection limit. However, since none of their additional five 
X-ray detected LBGs in the CDF-N (over and above the 
four LBGs detected in the CDF-N 2 Ms catalog by Alexan- 
der ct al. 2003 and discussed in Nandra et al. 2002) have 
luminosities L2-8 koV > 10^^ ergs s^^, this is not relevant 
to the issue of determining whether there is a decline in 
the L2-8 kev = lO'^^ - 1044-5 ergs s'^ range. 

Finally, N05 stated that their LBG color selec- 
tion mostly avoids any concerns about spectroscopic 
incompleteness, since the probability that the non- 
spectroscopically identified LBGs are z ^ 3 galaxies is 
about 96%. The high spectroscopic and photometric com- 
pleteness of our 2 = 2 — 3 X-ray sample mitigates this issue 
for our current analysis. 

Probably the most important concern about the N05 
methodology is that they must correct for two window 
functions (a very incomplete optical LBG selection func- 
tion that is not well understood for AGNs, and an X-ray 
selection hmction) rather than one. In order to determine 
whether their optical LBG selection function is the same 
for galaxies with substantial AGN contributions as it is for 



those without, they would need to spectroscopically iden- 
tify a complete X-ray sample. However, once they had 
undertaken such a spectroscopic survey to calibrate their 
optical LBG selection, then there would be no need for 
them to redo the AGN number densities using the LBG 
method. 

6. SUMMARY 

In summary, we constructed two uniform. X-ray flux- 
limited, highly spectroscopically complete subsamples of 
160 sources (eight of which are stars) with negligible X- 
ray incompleteness from the CDF-N 2 Ms X-ray data to 
compute AGN number densities in the z = 2 — 3 redshift 
interval for both high (L2-8 kcV = — 10^^) and in- 
termediate (L2-8 kcV = lO^^ - 1044) rest-frame 2-8 keV 
luminosity intervals. Of the 152 non-stellar sources, 102 
are spectroscopically identified. In this paper, we used new 
UH 2.2m ULBCAM J and H band imaging and Spitzer 
IRAC imaging to estimate IR photometric redshifts for 
the sources without spectroscopic redshifts, increasing our 
identified non-stellar sample to 144. Our final galaxy sub- 
samples contain no more than eight unidentified sources. 

We then calculated the z = 2 — 3 AGN number den- 
sities (for all spectral types together and for broad-line 
AGNs alone) for both luminosity intervals and compared 
them with those at z < 1.5, which we calculated from the 
spectroscopic sample of Barger et al. (2005). We find clear 
evidence for a decrease in the L2-8 kcV = lO^'^ — 10^^ AGN 
number densities at z > 1 and can reject the hypothesis 
that the number densities remain flat to z = 2 — 3 at above 
the 5a level. 
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